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ABSTRACT 

It is generally believed that the magnetic free energy accumulated in the corona serves as a 
main energy source for solar explosions such as coronal mass ejections (CMEs). In the frame- 
work of the flux rope catastrophe model for CMEs, the energy may be abruptly released ei- 
ther by an ideal magnetohydrodynamic (MHD) catastrophe, which belongs to a global magnetic 
topological instability of the system, or by a fast magnetic reconnection across preexisting or 
rapidly-developing electric current sheets. Both ways of magnetic energy release are thought to 
be important to CME dynamics. To disentangle their contributions, we construct a flux rope 
catastrophe model in the corona and solar wind and compare different cases in which we ei- 
ther prohibit or allow magnetic reconnection to take place across rapidly-growing current sheets 
during the eruption. It is demonstrated that CMEs, even fast ones, can be produced taking 
the ideal MHD catastrophe as the only process of magnetic energy release. Nevertheless, the 
eruptive speed can be significantly enhanced after magnetic reconnection sets in. In addition, 
a smooth transition from slow to fast eruptions is observed when increasing the strength of the 
background magnetic field, simply because in a stronger field there is more free magnetic energy 
at the catastrophic point available to be released during an eruption. This suggests that fast and 
slow CMEs may have an identical driving mechanism. 

Subject headings: Sun: corona — Sun: magnetic fields — Sun: coronal mass ejections(CMEs) 



1. INTRODUCTION 

It is generally believed that the accumulated 
magnetic free energy serves as a main energy 
source for the spectacular solar eruptive phenom- 
ena such as coronal mass ejections (CMEs) (see 
reviews by Forbes (2000) and Low (2001)), but 
it remains open how the magnetic energy is re- 
leased. Among various scenarios, the flux rope 
catastrophe mechanism is a very promising one 
(reviewed recently by Lin et al. (2003) and Hu 
(2005)). The catastrophe is an ideal magnetohy- 
drodynamic (MHD) process belonging to a global 
magnetic topological instability of the system. It 
releases energy without ohmic heating, especially 
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suitable for CMEs without associated flares (e.g., 
Forbes & Isenberg, 1991; Isenberg et al., 1993; 
Forbes & Priest, 1995; Hu et al., 2003b). A by- 
product of the catastrophe is the formation of one 
or more electric current sheets, which grows at the 
Alfvenic time scale and provides proper sites for 
fast magnetic reconnection. Such a reconnection 
further releases the stored magnetic energy and 
should be responsible for a solar flare associated 
with a CME event (e.g., Forbes & Lin, 2000; Lin 
& Forbes, 2000). Thus, in terms of this scenario, 
there are two means of magnetic energy release 
process involved in solar eruptions, which are the 
ideal MHD catastrophe and resistive magnetic re- 
connection. Both processes are thought to be im- 
portant to the CME initiation and acceleration. 
It is therefore necessary and important to disen- 
tangle their contributions to the CME dynamics 
and energetics. This article serves as a first step 
to solve this problem under the framework of the 
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specific catastrophe model for CMEs. 

The eruptive flux rope after catastrophe is ac- 
celerated mainly by the net or unbalanced Lorentz 
force, which is contributed by various source cur- 
rents both inside and outside the flux rope, in- 
cluding that along each electric current sheet and 
the corresponding image current, and the back- 
ground potential field which is determined by the 
normal component of the background field at the 
solar surface. To reveal the roles of these various 
components of magnetic forces in sustaining the 
flux rope in equilibrium and causing the eruption, 
Chen et al. (2006b) conducted detailed force bal- 
ance analyses. It was found that the primary lift- 
ing force is provided by the azimuthal current in- 
side the rope and its image below the photosphere, 
which is mainly balanced by the pulling force pro- 
duced by the background potential field when the 
rope is in equilibrium. During an eruption of the 
rope caused by the catastrophe, the force associ- 
ated with the rapidly-developing current sheet (s) 
constitutes an additional and significant restoring 
force that decelerates the rope. This force will be 
greatly reduced or even eliminated once magnetic 
reconnection sets in across the current sheet (s), 
as a consequence the magnetic reconnection may 
cause a further acceleration of the flux rope. This 
will be confirmed and quantified by the calcula- 
tions presented in this article with an axisymmet- 
ric flux rope catastrophe model. The effect of the 
solar wind plasma with an equatorial current sheet 
above the streamer cusp point extending to in- 
finity is considered so as to get a more realistic 
description of the CME acceleration and propaga- 
tion. Note that in most previous studies on the 
flux rope catastrophe this effect was not included. 
Sun & Hu (2005) did consider the effect of the solar 
wind on the rope catastrophe and the consequent 
eruption. They suggested that the flux rope catas- 
trophe can serve as a mechanism for slow CMEs. 
Yet, they investigated only the cases correspond- 
ing to a relatively weak background field in terms 
of ideal MHD. Based on their work, we take the 
effect of both the solar wind and magnetic recon- 
nection into account in the present study, which 
may be regarded as a starting point to disentan- 
gle the impacts of MHD catastrophe and magnetic 
reconnection on the CME dynamics in the solar 
wind. We put our focus on the differences between 
the flux rope dynamics in two situations in which 



we either prohibit or allow magnetic reconnection 
to take place across the rapidly-growing current 
sheets. We introduce the background corona and 
solar wind model and the flux rope catastrophe 
model in the following section. In § 3.1 we com- 
pare the solutions with and without magnetic re- 
connection to illustrate its impact on the CME 
dynamics. We emphasize the effect of the back- 
ground magnetic field strength in § 3.2 and provide 
a brief summary of this article with discussions in 
the last section of this article. 

2. FLUX ROPE CATASTROPHE MODEL 
IN THE CORONA AND SOLAR WIND 

Coronal magnetic flux ropes arc in essence 
three-dimensional structures with two ends an- 
chored to the photosphere. Considering that the 
length of a magnetic flux rope is much larger 
than its diameter, the axisymmetrical simplifi- 
cation (also called 2.5-D models) can be used to 
approximate the realistic flux rope system. With 
this approximation, the magnetic field B can be 
expressed with a magnetic flux function tp(t,r,6) 
in spherical coordinates (r, 6, <p), 

The derived resistive MHD equations in spheri- 
cal coordinates are written into the following form 
(see also Ding et al., 2006): 
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The symbols p, v, T, and j represent the density, 
flow velocity, temperature and the current density, 
respectively, p is the vacuum magnetic permeabil- 
ity, R the gas constant, G the gravitational con- 
stant, Mq the solar mass, and 7 the polytropic 
index which is taken to be 1.05 so as to obtain 
a reasonable solution of the background steady- 
state solar wind. The temperature and density at 
the coronal base are taken to be T = 2 x 10 6 
K and p = 1.67 x 1CT 13 kg m~ 3 , respectively. 
The magnetic flux function at the coronal base is 
taken to be the same as a dipolar field given by 
ip(r = Rq) = iposm 2 9/R Q where Vo = B R Q 2 . 
Bq represents half of the magnetic field strength 
at the polar hole (6 = 0) on the solar surface, 
which can be used to represent the strength of the 
background magnetic field of the system. Note 
that in this study the magnetic flux function at 
the solar surface is fixed to its initial distribu- 
tion which is determined by Bo. For the effect 
of the photospheric magnetic flux distribution on 
the coronal flux rope catastrophe, please check a 
recent paper written by Sun et al. (2007) . The pa- 
rameter B will be freely adjusted. For example, 
■00 = 9.7 x 10 13 Wb with B taken to be 2 G giving 
the case investigated by Sun & Hu (2005). It is 
defined that the ideal MHD situation corresponds 
to a zero magnetic resistivity ij. In the resistive 
situation, an anomalous homogeneous resistivity is 
used with r\ given by 77 = rj pv 8 oRQ where 770 = 0.1 
and v s0 = y/WJh = 181.8 km s" 1 . The steady- 
state polytropic solar wind solution is obtained by 
solving the above MHD equations with 77 = 0. The 
magnetic topology (white lines) and the velocity 
color contour map for Bq = 6 G from 1 Rq to 20 
Rq are illustrated in Figure la. For other values 
of B , the field topology and velocity distribution 
are basically similar to that shown in this figure. 
It can be seen that the solution is characterized by 
a typical streamer-current sheet-solar wind config- 
uration. The cusp point is located at about 3 Rq, 
and the flow velocity reaches up to 400 km s _1 at 
about 10 Rq along the equator. 

Based on the obtained corona and solar wind 
solution, we let a flux rope with prescribed mass, 
toroidal and poloidal magnetic fluxes and 5 , ip ) 
emerge from the equator at the coronal base. The 
detailed emerging process has been described pre- 



viously by Hu et al. (2003b) and Chen et al. 
(2006a) and will not be repeated here. Special 
numerical measures are taken to maintain *f/ p and 
^ v invariant and equal to their initial given values 
during the simulation (see Hu et al., 2003b). Fig- 
ure lb exemplifies the magnetic configuration and 
velocity contour map of a flux-rope system with 
the solar wind, where the border of the original 
flux rope is depicted with a green circle. It can be 
seen that the closed field region of the streamer ex- 
pands apparently with the emergence of the flux 
rope. Such a swelling of a coronal streamer is often 
observed with the white-light coronagraphs before 
CMEs (e.g., Howard et al., 1985). The mass con- 
tained by the flux rope per radian in the azimuthal 
direction is set to be ^-Mq where the unit of mass 
M = PoRq = 5.643 x 10 13 kg, the poloidal flux * p 
is taken to be 0.3 in units of tpo while the toroidal 
flux ^ v is changeable. Thus one may find MHD 
solutions with different values of ^ v to examine 
whether a catastrophe occurs, and find out the 
meta-stable state of the system characterized by 
^ p = 0.3 and a specific value of ^ v which de- 
pends on B - Starting from this state any slight 
increase of $> p or ^ v may excite the catastrophe. 
Therefore, the state is taken as the initial state 
for our simulation of the flux rope eruption. We 
choose to increase ^ p from 0.3 to 0.305 at t — so 
as to trigger the catastrophe. Physically speaking, 
the increase of the poloidal flux can be achieved 
by a twist of a long three-dimensional flux rope 
anchored to the photosphere. Besides the mass, 
magnetic fluxes of the flux rope, the helicity of the 
flux rope is also of interest to the study on CMEs. 
In an axisymmetric system like that investigated 
in this article, the two-dimensional magnetic he- 
licity can be calculated by the following integral 
according to Hu et al. (1997), 

Hj* = 27T J J tpB v rdrd6, (9) 

where the factor 2ir comes from the integral over 
the azimuthal direction, which should be removed 
if one wants to evaluate the magnetic helicity per 
radian in the azimuthal direction. 

When the flux rope breaks away from the sur- 
face and erupts upwards, a current sheet may de- 
velop below the flux rope. It is well known that 
numerical pseudoreconnection takes place across 
the current sheet in most numerical simulations, 
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which causes a false transfer of poloidal flux from 
the background to the flux rope and results in a 
topological change. In this work, we take special 
measure to prohibit such numerical reconnection 
in order to investigate the flux rope dynamics in 
the framework of ideal MHD. The magnetic flux 
function tp along the current sheet is invariant, 
which is known a priori, and any reconnection 
across the sheet reduces it in the present simu- 
lation. We therefore reassign ip along the current 
sheet to the known constant value at each time 
step. This technique, first proposed by Hu et al. 
(2003b), effectively eliminates numerical reconnec- 
tion across the equatorial current sheet. Note that 
this special measure is not employed in our calcu- 
lations in terms of resistive MHD with a non-zero 

T). 

The calculations are carried out in a domain 
of R Q < r < 30i? Q and < 9 < tt/2, which is 
discrctized into 150 x 90 grid points. The grid 
spacing increases according to a geometric series 
of a common ratio 1.024 along the radial direc- 
tion from 0.02 at the solar surface to 0.71 at the 
top boundary. And a uniform mesh is adopted in 
the 9 direction. The multistep implicit scheme de- 
veloped by Hu (1989) is used to solve the MHD 
equations. For the eruptive solutions, the calcula- 
tions are terminated once the top part of the ejecta 
reaches the upper boundary. 

3. Numerical results 

In this section, we first present and compare the 
solutions given by the ideal and resistive MHD cal- 
culations for the case with B = 6 G as a first step 
to disentangle the impacts of MHD catastrophe 
and magnetic reconnection on the CME dynamics 
in the solar wind. Then, we investigate the ef- 
fect of the background magnetic field strength by 
comparing results with different values of B - 

3.1. Impact of magnetic reconnection on 
flux rope dynamics 

As mentioned previously, our simulations on 
the flux rope catastrophe and eruption start from 
an equilibrium state which is a meta-stable flux 
rope system in the corona and solar wind back- 
ground. The catastrophe is triggered by a slight 
increase of the poloidal flux "Jp of the flux rope 
from 0.3 to 0.305 (or 8.7 - 8.85 x 10 13 Wb in phys- 



ical units) at t = with the critical axial flux in 
the rope * v = 0.209 (6.06 xlO 13 Wb) for B a =6 
G. After that, the flux rope starts to break away 
from the photosphere and erupts upwards. Figures 
lc and Id show the magnetic topology and veloc- 
ity color contours at the same instant (t = 280 
minutes) for the two solutions with and without 
magnetic reconnection (i.e., the cases with r\ ^ 
and 77 = 0). An apparent difference between the 
two solutions lies in whether a current sheet de- 
velops below the flux rope. The sheet forms and 
grows with the rope eruption in the ideal MHD 
case, while it is eroded by magnetic reconnection 
in the resistive MHD calculation. It is also ap- 
parent that a significant part of the poloidal flux 
has been transferred from the background to the 
flux rope, and a new streamer appears as the af- 
termath of magnetic reconnection in the resistive 
calculation. The rate of magnetic flux transfer is 
mainly determined by the effective resistivity con- 
sisting of the anomalous resistivity and the nu- 
merical resistivity involved in the pseudoreconnec- 
tion. Unfortunately, at this time, it is not pos- 
sible to eliminate the pseudoreconnection in our 
calculations in terms of resistive MHD. Therefore, 
it is difficult to control the flux transfer rate by 
simply adjusting the magnetic resistivity in this 
case. Further discussion regarding this issue will 
be given in our discussion section. Another major 
difference is the color distribution of the velocity 
contour maps which indicates how fast the flux 
rope ejecta is. It can be seen that the rope erupts 
faster in the solution with magnetic reconnection, 
as will be quantitatively revealed in Figure 2. 

In Figure 2, we plot the profiles of the helio- 
centric distance, velocity, and acceleration of dif- 
ferent parts of the flux rope ejecta, including the 
cusp point (in dotted), the rope top (in dashed), 
the rope axis (in solid), and the rope bottom (in 
dot-dashed) , left panels for the case without recon- 
nection and right for the reconnection case. It can 
be seen that in both cases the flux rope starts to 
take off at t f« 70 minutes. There is an apparent 
delay of about 50 minutes of the time when the 
cusp point starts to move upwards rapidly. The 
delay reflects the time taken for the eruptive flux 
rope to propagate from the coronal base to the ini- 
tial cusp point location in the corona. We can see 
that the cusp point undergoes the most dramatic 
acceleration in both cases: in about 70 minutes the 
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velocity of the cusp point reaches up to 800 km s _1 
with the maximum acceleration being 300 m s~ 2 
in the ideal case, and 1200 km s _1 and 500 m s~ 2 
for the resistive one. On the other hand, it takes 2 
to 3 hours for the flux rope to be accelerated to the 
maximum speed. After the maximum, the veloci- 
ties become more or less constant. The velocities 
of different parts of the ejecta vary significantly 
from the cusp point to the rope bottom. For exam- 
ple, at t — 280 minutes the exact moment at which 
the snapshot is taken for Figure 1, the velocity de- 
creases monotonically from about 800 km s _1 at 
the cusp point to 350 km s _1 at the rope bottom 
in the ideal case, and from 1100 km s _1 to 750 km 
s _1 in the resistive case. This monotonic decrease 
of velocities from the leading to trailing edges has 
been often observed by measurements of flux-rope 
like structures in the interplanetary space (e.g., 
Gosling et al., 1998), which simply indicates that 
the rope undergoes a rapid expansion during its 
eruption. The velocity profiles given by both so- 
lutions in Figure 2 as well as that in the following 
figure arc in a good agreement with a recent statis- 
tical study on the CME accelerations, which indi- 
cates that a CME usually undergoes multiphased 
kinematic evolution including an initial slow rise 
phase and a main rapid acceleration phase in the 
inner corona, and a relatively smooth propagation 
phase in the outer corona (Zhang & Dere, 2006). 
The magnitude and duration of the main accel- 
eration given by our calculations are also in line 
with their statistical results. Taking the resistive 
case shown in Figure 2 as an example, the main 
acceleration phase starts from t « 70 minutes and 
ends at about t = 190 minutes lasting for nearly 
2 hours. Comparing the solutions with and with- 
out magnetic reconnection, it is apparent that the 
speeds and accelerations of the flux rope are sig- 
nificantly enhanced in the case involved with mag- 
netic reconnection. Further discussions regarding 
the roles of magnetic reconnection in the CME dy- 
namics will be given in our discussion section. 

3.2. Effect of background field strength 

In this subsection, we present numerical results 
given by calculations with different values of the 
parameter Bq, which represents the strength of 
the background field and directly relates to the 
amounts of magnetic energy that can be stored 
and released in the system. In Figure 3, we plot 



radial profiles of the velocity and acceleration of 
different parts of the system including the cusp 
point (in dotted), and the rope top (in dashed), 
axis (in solid), and bottom (in dot-dashed) for the 
cases with (thick lines) and without (thin lines) 
reconnection, the left panels arc for the solution 
with Bq = 2 G and the right panels with Bq = 10 
G. Note that for clearness only the accelerations 
of the cusp point and rope axis are plotted in the 
lower panels. Similar to the above calculation the 
catastrophe is triggered by a slight increase of the 
dimcnsionless poloidal flux inside the rope ^ p from 
0.3 to 0.305 with the corresponding critical dimcn- 
sionless axial flux equal to 0.129 for the case 
with B = 2 G and 0.244 for B = 10 G. Note 
that the magnetic fluxes in physical units are listed 
in Table 1. Since the relative magnitude of the 
rope poloidal flux does not change with varying 
Bq, the size of the flux rope does not change ap- 
parently either. To be quantitatively, we checked 
the heliocentric distance of the rope axis v a which 
can be used to represent the size of the flux rope. 
It was found that r a = 1.42, 1.52, and 1.54 solar 
radii for the cases with B =2, 6, and 10 G, respec- 
tively. As a result, the plasma density inside the 
flux rope gets slightly different in different cases 
since the total mass contained by the rope does not 
change with Bq. In Table 1, we also list the two- 
dimensional magnetic helicity per radian in the 
azimuthal direction Ht calculated with Equation 
(9). Similar as the rope fluxes, Ht also increases 
dramatically with increasing Bq. We see that the 
physical properties of the flux rope differ signifi- 
cantly from case to case in our calculations, which 
naturally have impacts on the flux rope dynam- 
ics. Further discussion along this direction will be 
presented at the end of the following paragraph. 

It can be seen from Figure 3 that the most ob- 
vious difference between the two sets of solution 
is the magnitude of velocity and acceleration at 
different points of the ejecta. In the resistive cal- 
culation with Bq = 2 G, the main acceleration 
phase of the flux rope lasts for about 4 hours from 
t « 100 minutes to t « 350 minutes with speeds 
rising up to 670 km s _1 at the cusp point and to 
460 km s _1 at the rope bottom, which is followed 
by the so-called propagation phase with a nearly 
constant speed (Zhang & Dere, 2006). In the ideal 
calculation with the same value of Bq, the veloc- 
ity keeps increasing till t — 600 minutes, yet the 
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corresponding acceleration gets smaller than 5 m 
s~ 2 . In the resistive case with Bo = 10 G, the 
main acceleration phase lasts for about 2 hours 
starting from t w 40 minutes and ending at about 
t = 150 minutes with the maximum acceleration 
rising up to 750 m s~ 2 . After the maximum of 
velocity is reached, which is about 1600 km s _1 
for the cusp point and 1150 km s _1 for the rope 
bottom, the flux rope gets decelerated gradually 
to a velocity of 1300 km s _1 at the cusp point and 
900 km s _1 at the rope bottom. In the ideal case 
with B = 10 G, the main acceleration phase also 
lasts for about 2 hours with the maximum accel- 
eration being about 430 m s -2 for the cusp point 
and 220 m s~ 2 for the rope axis. In the propa- 
gation phase following the velocity maximum, the 
velocities decrease slightly. It can be seen that 
CMEs, even fast ones, can be produced taking the 
ideal MHD catastrophe as the only process of mag- 
netic energy release. It is also true, again, that 
the eruptive speeds are significantly enhanced af- 
ter magnetic reconnection sets in. We point it out 
in passing that a smooth transition of eruptions 
from slow to fast can be obtained when varying 
Bo continuously with a stronger background field 
corresponding to a faster eruption, in line with the 
very recent study on the effect of photospheric flux 
distribution on the flux rope dynamics by Sun et 
al. (2007). The physical cause of such a behav- 
ior can be easily understood from the following 
simple energy analysis. Since the pattern of the 
magnetic flux distribution at the coronal base re- 
mains the same for all cases we have discussed, 
the associated open field energy must be propor- 
tional to the square of B . On the other hand, 
the percentage by which the catastrophic energy 
threshold exceeds the open field energy varies in a 
much smaller range. It reads 6.1%, 8.8% and 9.6% 
for Bo = 2 G, 6 G and 10 G, respectively. Such a 
result is consistent with previous similar calcula- 
tions (e.g., Hu et al., 2003b, Li & Hu, 2003; Chen 
ct al., 2006a). Therefore, the total ammount of 
magnetic free energy of the system at the catas- 
trophic point is mainly determined by the overall 
strength of the background field in spite of the dra- 
matic differences in the flux rope properties at the 
catastrophic point as listed in Table 1. In sum- 
mary, we argue that the stronger the background 
field is, the more magnetic free energy is available 
for the flux rope system at the catastrophic point, 



which leads to a faster eruption of the flux rope. 

To shed more light on the effect of magnetic 
reconnection on the CME dynamics, we calculate 
the total increase in kinetic energy of the system 
compared with that of the pre-eruption state (i.e., 
the state at t = 0), represented by AEk- Figure 
4 shows the temporal profiles of AEk per radian 
in the azimuthal direction in units of 5.38 x 10 31 
ergs for the three sets of solutions with B =2 G 
(in solid), 6 G (in dotted), and 10 G (in dashed), 
where thick and thin lines represent the results 
with and without magnetic reconnection. The ve- 
locity and acceleration for these solutions have 
been illustrated in Figures 2 and 3. We can see 
that AEk tends to reach an asymptotic value in 
all cases. For each set of solutions, the asymptotic 
value of AEk in the resistive case is about 2 to 3 
times of that in the ideal case indicating that the 
MHD catastrophe and magnetic reconnection, the 
two means of magnetic energy release process, are 
of comparable importance on the CME accelera- 
tion for the resistive MHD situation. 

4. Conclusions and discussion 

In terms of the catastrophe theory, there are 
two main processes energizing the solar erup- 
tions: MHD catastrophe and magnetic reconnec- 
tion. This article serves as a first step to disen- 
tangle their contributions to the CME dynamics. 
To do this, we construct a flux rope catastrophe 
model in the corona and solar wind and compare 
different cases in which we either prohibit or al- 
low magnetic reconnection to take place across 
rapidly-growing current sheets during the erup- 
tion. For simplicity, a polytropic process with the 
polytropic index 7 =1.05 is used to produce the 
background solar wind solution. The catastrophe 
and the consequent eruption is triggered by a tiny 
increase of the rope poloidal flux, which reflects a 
slight twist of the ends of a long three-dimensional 
realistic flux rope anchored to the photosphere. It 
is demonstrated that CMEs, even fast ones, can 
be produced taking the ideal MHD catastrophe 
as the only process of magnetic energy release. 
Nevertheless, the eruptive speed can be signifi- 
cantly enhanced after magnetic reconnection sets 
in. In addition, a smooth transition from slow 
to fast eruptions is yielded when increasing the 
strength of the background magnetic field, i.e., a 
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stronger field, in which more free magnetic energy 
gets available at the catastrophic point, enables a 
faster eruption. This suggests that fast and slow 
CMEs may have an identical driving mechanism. 

Based on previous and present studies taking 
catastrophe as the principle driving mechanism 
of CMEs, we argue that the MHD catastrophe 
is probably the main means of energy release for 
CMEs at least in the initial phase. It releases en- 
ergy without ohmic heating and provide accelera- 
tions with the Lorentz force, especially suitable for 
non-flare associated CMEs. A by-product of the 
catastrophe is the formation of one or more elec- 
tric current sheets, which proceeds at the Alfvenic 
time scale and produces conditions favoring fast 
magnetic reconnection. Such a reconnection, if 
takes place, further releases the magnetic energy 
through the following two aspects. Firstly, the 
magnetic energy is converted into thermal and ki- 
netic energies of plasma particles at the reconnec- 
tion site. This process is believed to account for a 
solar flare associated with a CME event. Secondly, 
the restoring force contributed by the current in 
the current sheet is significantly reduced or even 
eliminated, and the magnetic topology changes 
with the magnetic reconnection. This also pro- 
duces a significant acceleration of the flux rope in 
addition to that caused by catastrophe. It can 
be seen from our quantitative calculation that the 
MHD catastrophe and magnetic reconnection, the 
two magnetic energy release processes, may have 
comparable impacts on the CME dynamics during 
the main acceleration phase of CMEs. 

A major subject of this article is to estimate 
the impact of magnetic reconnection on the flux 
rope dynamics. For this purpose, we compared 
solutions given by calculations with and without 
magnetic reconnection. However, this work suf- 
fers from the facts that the anomalous magnetic 
resistivity is artificially given and the numerical 
pseudoreconnection is unavoidable in the resistive 
calculations. Since the effective resistivity includ- 
ing both the anomalous and the numerical ones 
is believed to be much larger than the realistic 
value in the corona and solar wind, the velocity 
profiles given by our study for the case with mag- 
netic reconnection should be taken as the upper 
bound for the realistic situation. As mentioned 
in the text, the effective resistivity is a crucial 
factor determining the transfer rate of magnetic 



flux from the background to the flux rope. Obser- 
vationally, this transfer rate can be evaluated by 
extrapolating the photospheric field to the corona 
and counting the change of the total magnetic flux 
in the aftermath of a CME event, e.g., in the coro- 
nal dimming region (see, e.g., Jing et al., 2005 and 
Qiu & Yurchyshyn, 2005). Further theoretical en- 
deavor should utilize these relevant observational 
constraints on the flux transfer rate in the mod- 
elling of a specific event. 

Our study on the effect of magnetic field 
strength reveals a smooth transition from slow 
to fast eruptions when increasing the background 
field strength. This is in support of the argument 
that slow and fast CMEs may be driven by a single 
identical mechanism, also in line with recent sta- 
tistical results contradicting with the traditional 
bimodal classification of slow (gradual) and fast 
(impulsive) CMEs (e.g., Sheeley et al., 1999; An- 
drews & Howard, 2001). For instance, it is shown 
by Yurchyshyn et al. (2005) and Zhang & Dere 
(2006) that the velocity and acceleration of a large 
amount of CME events have a continuous distribu- 
tion instead of a bimodal one, and by Vrsnak et al. 
(2005) that flare- and nonflare- associated CMEs 
have quite similar characteristics in the LASCO 
C2 and C3 fields of view. The velocity profiles for 
slow CMEs given by our model, say, the solution 
corresponding to a weak background field, show 
a gradual acceleration, while that for fast CMEs 
present a rapid acceleration and a discernable de- 
celeration following the main acceleration phase. 
This behavior is probably a result of the coupling 
process with the background solar wind plasma 
according to our preliminary analysis. When the 
magnetic energy released in an eruption is not 
enough to accelerate the flux rope ejecta to the 
speed of the background plasmas, the ejecta may 
get gradually accelerated and gain more energy 
through the coupling to the solar wind. On the 
other hand, the ejecta may get decelerated and 
lose energy through similar coupling process with 
the solar wind while the released magnetic energy 
is enough to push the flux rope outwards with a 
velocity faster than that of the background. 

There exist contradicting discrepancies between 
the present thick-rope model in axisymmetrical 
spherical geometry and that published in the liter- 
ature in terms of thin-rope models in 2-D cartesian 
geometry. Firstly, an infinite amount of energy is 
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required to open up a closed magnetic field in 2-D 
Cartesian geometry (Hu et al., 2003b), therefore, 
it is energetically impossible to open the overly- 
ing field and to let the flux rope escape to infinity 
without magnetic reconnection, as demonstrated 
by the catastrophe models assuming Cartesian ge- 
ometry (e.g., Lin & Forbes, 2000). On the other 
hand, in the spherical geometry the corresponding 
open-field energy is finite and it can be exceeded 
by the flux rope system as already shown by many 
calculations (e.g., Choe & Cheng, 2002; Hu et al., 
2003b; Li & Hu, 2003; Flyer et al., 2004; Sun & 
Hu, 2005; Zhang et al., 2005; Peng & Hu, 2005; 
Ding & Hu, 2006; Chen et al., 2006a). Thus, mag- 
netic reconnection may not be necessarily required 
for the flux rope to get escaped from the Sun in the 
spherical model. Secondly, as pointed out previ- 
ously, the self-interaction of the azimuthal current 
inside the flux rope by itself results in an outward 
radial force on the rope, which comes from the 
curvature of the rope surrounding the Sun (Chen, 
f989; Lin et al., 1998; Krall et al., 2000; Chen 
et al., 2006b). This self-force, together with that 
contributed by its image current below the pho- 
tosphere, serves as a dominant driving force for 
the rope eruption. However, in the 2-D Carte- 
sian models this self-force is trivially zero by the 
symmetry of an infinitely long straight current, 
this gives another basic difference between the 2- 
D Cartesian and spherical models. 

No matter what geometry is used, so far most 
flux rope models have been limited to 2-D anal- 
yses, as a necessary simplification for practical 
tractability. Yet, the two ends of a 3-D flux rope 
are believed to be anchored to the solar surface, 
ft remains open regarding how the catastrophic 
behavior of the flux rope may change under this 
situation. Finally, since how the corona and so- 
lar wind plasma is heated and accelerated is still 
a big issue to be resolved, the polytropic process 
is assumed conveniently to obtain the reasonable 
background solar wind solution in this work. How- 
ever, it should be noted that the polytropic solar 
wind solution is too simple to account for some re- 
alistic properties of the solar wind. For example, 
the effect of the fast solar wind is not included in 
this study. This will certainly affect the propaga- 
tion of the ejecta at high latitudes, yet may not 
be very important to the study on the CME prop- 
agation along the equatorial plane. In future we 



consider to employ sophisticated heating functions 
to produce a more realistic solar wind background 
(see, e.g., Chen &Hu, 2001; Hu et al., 2003a) for 
a more elaborated study on the CME propagation 
in the meridional plane. 

This work was supported by grants NNSFC 
40404013, NSBRSF G2006CB806304, and NNSFC 
10233050 in China. 
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(Wb) 2.91 xlO 13 8.70 xlO 13 1.45 xlO 
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H T (Wb 2 ) 1.45 xlO 27 2.12 xlO 28 6.32 x10 
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This 2-column preprint was prepared with the AAS IATjrjX 
macros v5.2. 
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Fig. 1. — The magnetic topology (white lines) and 
velocity color contour map for Bq — 6 G from 1 Rq 
to 20 Rq for: (a) the background corona and so- 
lar wind solution before the emergence of the flux 
rope, (b) the pre-eruption state of the flux-rope 
system in the solar wind background, (c) the solu- 
tion with an erupting flux rope at t = 280 minutes 
without magnetic reconnection, and (d) the erup- 
tive solution with magnetic reconnection at the 
same instant as panel (c). The outer boundary 
of the original flux rope is depicted with a green 
circle. 



Fig. 3. — The radial profiles of the velocity and 
acceleration of different parts of the system in- 
cluding the cusp point (in dotted), and the rope 
top (in dashed), axis (in solid), and bottom (in 
dot-dashed), the left panels are for the case with 
B = 2 G and the right panels with B — 10 G 
for the cases with (thick lines) and without (thin 
lines) reconnection. 
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Fig. 2. — The temporal profiles of the heliocen- 
tric distance, velocity, and acceleration of differ- 
ent parts of the flux rope system, including the 
cusp point (in dotted), the rope top (in dashed), 
the rope axis (in solid), and the rope bottom (in 
dot-dashed). Left panels are for the case without 
reconnection and right for the reconnection case. 



5 ; 

/ 

4 i I 




100 200 300 400 500 600 
Time (minutes) 



Fig. 4. — The temporal profiles of the total in- 
crease in kinetic energy of the system AEk in units 
of 5.38 x 10 31 ergs (per radian in the azimuthal 
direction) for three sets of solutions (Bo— 2 G (in 
solid), 6 G (in dotted), and 10 G (in dashed)) with 
(thick lines) and without (thin lines) magnetic re- 
connection. 
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